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Abstract — Interference alignment is a transmission teclinique 
for exploiting all available degrees of freedom in the frequency- 
or time-selective interference channel with an arbitrary number 
of users. Most prior work on interference alignment, however, 
neglects interference from other nodes in the network not 
participating in the alignment operation. This paper proposes 
three generalizations of interference alignment for the multiple- 
antenna interference channel with multiple users that account for 
colored noise, which models uncoordinated interference. First, 
a minimum interference-plus-noise leakage (INL) algorithm is 
presented, and shown to be equivalent to previous subspace meth- 
ods when noise is spatially white or negligible. This algorithm 
results in orthonormal precoders that are desirable for practical 
implementation with limited feedback. A joint minimum mean 
squared error design is then proposed that jointly optimizes the 
transmit precoders and receive spatial filters, whereas previous 
designs neglect the receive spatial filter. Finally, a maximum 
signal-to-interference-plus-noise ratio (SINR) algorithm is de- 
veloped that is proven to converge, unlike previous maximum 
SINR algorithms. The sum throughput of these algorithms is 
simulated in the context of a network with uncoordinated co- 
channel interferers not participating in the alignment protocol. 
It is found that a network with co-channel interference can benefit 
from employing precoders designed to consider that interference, 
but in extreme cases, such as when only one receiver has a large 
amount of interference, ignoring the co-channel interference is 
advantageous. 



I. Introduction 

Interference channels model a network of simultaneously 
communicating node pairs. In these channels, each transmitter 
has data to send to only one receiver, which also observes 
interference from the other transmitters in the network. Anal- 
ysis of interference channels has shown that interference is 
not a fundamental limitation. In particular, with any sized 
interference channel with any number of users, the capacity for 
any given user will scale at half the rate of its interference-free 
capacity in the high transmit power regime 

The key to achieving a linear capacity scaling is interference 
alignment (lA) f2\, f^. With lA, interfering transmitters pre- 
code their signals to align in the unwanted users' receive space, 
allowing these receivers to completely cancel more interferers 
than they otherwise could. The signals can be aligned in any 
dimension, including time, frequency, or space. This can be 
viewed as a cooperative approach because the transmitters 

The authors are with the Department of Electrical and Computer Engi- 
neering, 1 University Station C0806, University of Texas at Austin, Austin, 
TX, 78712-0240 (email: {speters,rheath}@mail. utexas.edu, phone: (512) 471- 
1190, fax: (512) 471-6512). 

This work was supported by the DARPA IT-MANET program. Grant 
W91 lNF-07-1-0028, and a gift from Huawei Technologies, Inc. 

Copyright (c)2010 IEEE. Personal use of this material is permitted. How- 
ever, permission to use this material for any other purposes must be obtained 
from the IEEE by sending a request to pubs-permissions@ieee.org. 



neglect the performance of their own link to allow other 
users to perfectly cancel interference. This is in contrast to 
a provably suboptimal "selfish" approach where a transmitter 
ignores the interference it causes and aims simply to maximize 
its own data rate |4|. Interference alignment has been shown 
to achieve the maximum capacity scaling, also known as 
degrees of freedom, of the 7^-user interference channel, but 
at finite transmit power it offers suboptimal achievable sum 
rate. Consequently, there is interest in finding precoders for the 
interference channel that relax the perfect alignment constraint 
with the objective of obtaining better nonasymptotic sum rate 
performance. 

Alternative lA precoder designs have been proposed for 
the single-antenna interference channel with time or frequency 
selectivity [5], |6|. Closed-form lA precoders and achievable 
degrees of freedom for the multiple-input multiple-output 
(MIMO) interference channel with infinitely selective channels 
have also been found for some asymmetric antenna arrange- 
ments I?). Interference channels where precoding can only be 
done over one transmission slot are said to have constant or 
static coefficients. In this case, the degrees of freedom with 
linear precoding are unknown but have been hypothesized to 
be less than that with infinite selectivity |9|, while non- 
linear precoding might achieve KAI/2 degrees of freedom 
with M antennas at each transmitter and receiver fTol , | fTT| . 

A challenge in constant coefficient MIMO interference 
channels is that closed form solutions have been found in 
only a few special cases fl]. Algorithmic techniques, such 
as alternating minimization [121, have been proposed to find 
precoders and explore possible degrees of freedom for the 
general case p3|-p3|. Such algorithms are promising both for 
their ability to provide precoder solutions in a practical setting 
and their flexibility in application to arbitrary networks for 
which closed-form solutions are unknown. The subspace al- 
gorithms of |T3^-fT5^, however, still use alignment as the main 
objective, which is asymptotically optimal for the interference 
channel but has suboptimal throughput at finite SNR and other 
regimes. They also neglect colored noise, possibly caused by 
co-channel interference from outside the coordinating nodes. 
A maximum-SINR algorithm was proposed in [12], but this 
algorithm does not optimize a global objective, assumes white 
Gaussian noise, and is not shown to converge. 

In this paper we propose several alternative linear precoding 
designs for MIMO interference channels. While maximizing 
the sum rate is the primary objective, we do not directly 
maximize sum rate due to analytical intractability. Instead 
we approximate a sum rate maximization via algorithms with 
varying performance and complexity tradeoffs. First, we derive 
a generalization of subspace alignment that includes colored 
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noise, which biases the preferred alignment subspaces. The 
resulting objective, which minimizes the interference plus 
noise leakage (INL), results in orthogonal precoders amenable 
to quantized CSI. This algorithm is shown to be a special 
type of minimum mean squared error (MMSE) design, and 
at high SNR or white noise at all receivers is shown to 
reduce to the lA subspace methods of |jT3), (H). From this, 
interference alignment is shown to be an MMSE-type solution 
at infinite SNR, where interference-suppression filters are 
optimal. As with previous forms of interference alignment, the 
proposed minimum INL algorithm does not consider the signal 
power at any given user and is thus suboptimal with finite 
transmit power Further, this algorithm and previous designs 
derive precoders but neglect receiver design, which could be 
optimized jointly with the precoders. 

Inspired by the connection between mutual information and 
mean-square error fW\, we derive an explicit joint MMSE pre- 
coder/receiver design for the interference channel. Although 
it does not directly maximize the sum rate, the joint-MMSE 
design results in a higher sum rate than subspace methods. 
It does not lead to orthonormal precoders, making quantized 
feedback design more difficult. The MMSE design is shown to 
be a generalization of previous approaches for point-to-point 
and multiuser settings |17J-|21J. Further, the design is more 
computationally complex and requires more iterations at high 
SNR than subspace designs. MMSE-based designs have also 
recently been developed independently in p2) , ||23 |. 

To more directly optimize the sum rate we formulate a 
maximum SINR algorithm, which is proven to converge via 
alternating minimization of a global performance function. The 
maximum SINR algorithm derived in |jT3] is shown to be an 
approximation to the one derived in this paper. On average, the 
two are shown to have the same performance, but for any given 
channel realization may result in different sum rates. This 
design often has increased throughput relative to MMSE and 
subspace approaches, but finds nonorthogonal precoders and 
requires more channel state information if run in a distributed 
manner. 

In summary, this paper proposes three algorithms that span 
the tradeoff between performance and complexity for the static 
MIMO interference channel. The minimum INL algorithm 
has the same complexity as previous work but has improved 
performance when colored noise exists at any receiver The 
joint-MMSE design has further rate enhancements regardless 
of the noise covariance matrices but has a computationally 
more complex optimization and non-orthogonal precoders. 
The maximum SINR design has the best overall performance 
of all proposed strategies (in most cases, as shown in the 
simulations), but requires more channel state information than 
the previous designs and also results in non-orthogonal pre- 
coders which are difficult to quantize in a practical setting p4) . 
The proposed algorithms are then simulated alongside existing 
methods in regimes previously unconsidered in the literature. 
For example, the algorithms are simulated in an environment 
with an uncoordinated interferer that is not participating in 
the alignment protocol. This colors the noise at each receiver, 
and if its power is scaled with the rest of the transmitters, 
resulting in reduced capacity scaling. Each of the algorithms 



can outperform the others in different regimes, and each of 
these regimes is simulated and enumerated. 

The rest of this paper is organized as follows: Section III 
presents the system model under consideration; Section IlII 
presents the new MMSE and INL algorithms and derives 
a maximum SINR algorithm with proven convergence and 



analyzes each of the methods; Section IV presents simula- 
tions under uncoordinated interference and colored noise; and 
Section [V] concludes the paper and gives directions for future 
work. 

Before proceeding, we introduce notation. The log refers 
to log2. Bold uppercase letters, such as A, denote matrices, 
bold lowercase letters, such as a, denote column vectors, and 
normal letters a denote scalars. The letter E denotes expecta- 
tion, C is the complex field, M{a} is the real component of 
complex scalar a, min{a, b} denotes the minimum of a and b, 
i^^iij (A) is the matrix whose columns are the eigenvectors 
corresponding to the R smallest eigenvalues of matrix A, 
tr (A) is the trace of matrix A, \a\ is the magnitude of the 
complex number a, ||a|| is the Euclidean norm of vector a, 
and |A| is the determinant of square matrix A. A* is the 
Hermitian transpose of matrix A and A^^ is its inverse. The 
matrices I and are the identity matrix and all zero matrix, 
respectively, of appropriate dimension. Finally, we use {F^} 
when referring to the set of precoders and when referring to 
the precoder at transmitter £, and similarly for receive spatial 
filters {Gk}. 

II. System Model 

Consider the K-usei MIMO interference channel illustrated 
in Figure [T] with K transmit-receive pairs. A wireless channel 
links each receiver to each transmitter, but a given transmitter 
only intends to have its signal decoded by a single receiver 
The kth transmitter possesses antennas with which to 
transmit Sk < spatial streams, and the kth receiver (which 
is to decode the signal from the A:th transmitter) possesses 
Nk > Sk antennas. In some analysis and simulations, all 
users will have the same antenna configurations so that Mk — 
M, Nk = N, and Sk = S, Vfc; we denote this symmetric case 
as an (Af , N, K) interference channel with S streams per user 

This paper considers the narrowband MIMO interference 
channel where each link is static for the duration of a transmis- 
sion, but may change between successive transmissions. This 
is the block fading model where all the links in the network 
are constant for the period of transmission, creating a tractable 
approximation to more realistic continuous fading models. 
Linear precoding is done independently over each channel 
realization, favoring simplicity over the possible degrees of 
freedom gained by jointly precoding over realizations. This 
is the same model as previous work on algorithms for the 
interference channel |jT3)-|(l5). The transmission of all K 
users is synchronized such that each begins and ends each 
transmission simultaneously, and no frequency offsets exist in 
the network. We therefore take the standard approach p3) , 
1 15 1 and focus on the transmission of a single vector symbol 
Sfe from transmitter k G {1,...,K}, neglecting any time 
dependency. 
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where 



Interfering Link 
^ Communicating Link 



Fig. 1. The MIMO interference cliannel. Eacli transmitter is paired with a 
single receiver, and all links are non-negligible. 



Transmitter k uses linear precoder e (^A/^xSfc jq j^^^p 
Sk symbols in to its Mk transmit antennas. 



(1) 



where the transmitted symbols are i.i.d. such that Es^s^ — I, 
the precoder is normalized such that ||Ffc|||, < p^, and is 
the transmit power at transmitter k. Receiver k observes the 
signal 



K 



(2) 



£=1 



where ^ is the channel between transmitter i and receiver k 
and Vfc is Gaussian noise at receiver k with spatial covariance 
matrix Rj. = Evj-v^. For the analysis in this paper we assume 
that the channels {Hfc £} are each full rank and mutually 
independent, the transmitters send independent data (Es^s^ — 
for /c 7^ £) and all transmitted signals are statistically 
independent from the noise at any receiver (Es^v^. = for all 
{k, £) G {1, . . . , K}^). No assumptions are made on the noise 
power or covariance at any receiver. Rewriting Q, receiver k 
sees 



K 



Yk 



— Hfe^fcFfcSfe + ^ Hfc^^F^Sf + Vfc. 



(3) 



i=i 

l^k 



The vector s^. is the signal to be decoded by receiver fc, and the 
summation term in ([3]) is called coordinated interference, since 
it is caused by transmitters that may coordinate to minimize 
its effect. Once the precoders are designed, the instantaneous 
sum rate of the system is 



fc=i 



iTT T7I TP*TT* 



(4) 



K 



(5) 



£=1 



is the interference plus noise covariance matrix at receiver 
k. The instantaneous sum rate is an important metric for 
multiuser systems because of its ability to capture the total 
network throughput in a single scalar. Notice that i?sum 
assumes ideal non-linear decoding of the signal. Although 
the proposed algorithms of Section [III] will design linear 
processing matrices that can form a part of a linear receiver 
design, they serve mainly to simplify the optimization and 
design of the precoders. Design of high performance linear 
receivers is left to future work, except for the MMSE design 



presented in Section III-C Thus for fair comparison, the 
sum rate equations assume an ideal decoding for all precoder 
designs. 

Previous authors have shown that KM/2 spatial degrees of 
freedom are achievable in an (M, M, K) interference channel 
Degrees of freedom d is 

Csum (SNR) 



d = lim 

SNR-i-c 



log SNR 



(6) 



where Csum is the sum capacity of the network, rather than 
the sum rate for our linear precoding model presented in Q. 

The key idea of interference alignment is to make X^^^^fc 
interferers appear as Nk — Sk interferers at receiver k for each 
k by having them span a subspace of dimension — Sk of 
the A'fc -dimensional receive space. Mathematically, 



E 

l^k 



Hfc,£F. 



^ a,c[.*\VA;, 

i=l 



(7) 



where {c[*''} are basis vectors for the subspace at receiver 
k in which all interference must lie. Then receiver k can 
then resolve its Sk streams with a linear receiver interference- 
free [|T]. For the three user interference channel it is possible to 
directly find closed-form solutions to {F^ }. for any S < M/2. 
Such solutions for obtaining KM/2 degrees of freedom, 
however, in the {M, M, K) interference channel with K > Z 
users do not appear to be possible. Closed-form solutions, even 
for a reduced multiplexing gain, are unknown I^SJ except in 
special cases A viable alternative for the general case 
are alternating minimizations. The next section reviews the 
existing designs and proposes new algorithms for finding high- 
rate solutions at finite-SNR in the MIMO interference channel. 



III. Iterative Algorithms Via Alternating 
Minimization 

This section presents iterative solutions for precoders in the 
MIMO interference channel using an alternating minimization 
to solve various optimization objectives. This section proposes 
three new metrics which aim to approximate a sum rate 
maximization with better finite-SNR rates than previous work. 

The algorithms presented may be implemented in a dis- 
tributed or centralized manner similar to p3| . These algo- 
rithms share a common structure. Each algorithm is designed 
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to optimize a global objective J that incorporates the perfor- 
mance of each data link in the network. The objective is a 
function of the precoders {F^}, the channels {Hfc^} between 
all node^ and a processing matrix at each receiver, the 
structure of which will vary across design^ The free variables 
are the K precoders and K receive processing matrices. 

A closed-form solution for a global optimization of any 
of the objectives in this section is unknown. We therefore 
turn to an alternating minimizatiorj^ approach for the 2K 
variables |12|. In general, an alternating minimization ar- 
bitrarily initializes 2K — 1 variables and, assuming these 
variables are fixed, solves for the remaining one. It stores 
this solution, and moves to another variable, finding a new 
solution for it assuming the rest of the variables are fixed. 
Each variable in turn is solved for during each iteration. Note 
that this procedure is convenient only if there is a simple or 
even closed-form solution for each of the variables assuming 
the rest are fixed. Finally, for each of the designs, with 
the exception of the proposed maximum SINR design, the 
precoders may be derived in parallel, since their solutions at 
any step of the algorithm do not depend on each other. 



Although we could enforce a total power constraint on the 
precoders (and, coincidentally in this case, get the same solu- 
tion), orthogonality is desired in MIMO precoding designs to 



aid with feedback of channel state |25|. The receive subspace 
bases {^fc} are orthonormal by definition so that ^^^fe = I. 
The objective is thus 



minimize Jia({F£},{*}) 
subject to F|F£ = {l,...,is:} 



(10) 



*fe*fe 



I,fce{l,...,if}. 



The optimization ( 10 1 is intuitively pleasing since, with perfect 
interference alignment, JTia = 0, and without interference 
alignment, JTia > 0. That is, interference alignment, if 
possible, achieves the global minimum for this function. 
Deriving a closed-form solution to (lOi for if > 3 users 



is difficult due to the inter-dependence of each precoder 
and receive interference-free subspace. A simple approach, 
which is guaranteed to converge, is to use an alternating 
minimization y_2J. The derivation of this solution is in |13| 
and our previous work in p5| and is not included here for 
efficiency. At each step, the solution for each F^ is 



A. Subspace Optimization 

A direct algorithm for the interference channel inspired 
by interference alignment |rj is to precode the signal at 
transmitter i such that the coordinated interference caused 
by transmitter ^ at receiver fc 7^ ^ is nearly orthogonal 
to a subspace (with orthonormal basis $/j) of its receive 
space p3| , p3] . This subspace is then jointly designed along 
with the precoders to optimize an appropriate cost function. 
One way of performing this optimization is to minimize the 



K 



total "leakage interference" 1 13| that remains at each receiver 



after attempting to cancel the coordinated interference by left- 
multiplication with for each fc. The global function to 
optimize is thus 



JiA 



K 

E- 

fe=i 



E 



K 



*^^Hfc,,F,s, 



(8) 



The expectation in Jia and all subsequent analysis is over 
(and Vfe where applicable), k g {l,...,if}. Evaluating the 
expectation and exploiting independence of the signals. 



K K 



II, 



(9) 



k=l 1=1 



which is termed "interference leakage" in fT3^. The precoders 
{F^} are constrained to have mutually orthogonal columns 
with a per-stream power constraint so that F^F^ = ^I, V^. 

'Some of the algorithms will not make use of the data links, instead 
focusing on minimizing post-processing interference. 

^These matrices are not necessaiily designed to function as spatial equal- 
izers, instead serving mainly to simplify the design and optimization of the 
precoders. With the exception of our joint-MMSE algorithm, design of high 
performance linear receivers is left to future work. 

^Alternating maximizations can be converted into alternating minimiza- 
tions, so we focus on alternating minimization. 



(11) 



k=l 

\k^e 



and, with all precoders given, the solution for each is 
(k 

^Hfc,,F,F|H^,, 1 . (12) 

\i^k 



To run the algorithm, arbitrary receive subspaces for each 
receiver are used for initialization and an arbitrary orthonormal 
basis for each subspace is found. This subspace is ideally 
reserved for user k's signal, thus coordinated interference 
at receiver k is ideally orthogonal to this subspace. Then, 
for each i, the algorithm finds the precoder matrix F^ such 
that total coordinated interference caused at each node (other 
than at node £) has maximum squared Euclidean distance 
between it and the subspace spanned by the columns of each 
using (111. Given these new precoders, the algorithm 
can update the receive subspaces to be those that span the 
columns of the matrices with minimum sum squared Euclidean 
distance to the interference caused by the fixed precoders 
using ( [T2[ i [15]. This can be carried out until JiA{t) < e 
if feasibility conditions are met, or JiA{t — 1) — JiA{t) < e 
otherwise, for an arbitrary convergence threshold e. 

Note that each receiver must still separate the desired spatial 
streams after the coordinated interference has been canceled 
with left multiplication of Standard linear designs, such 
as zero forcing or MMSE, can be employed for this purpose. 
Thus, the receiver can form a linear receive filter by 
multiplying ^k and the linear spatial filter W^, which ne- 
glects coordinated inter-user interference and equalizes only 
the desired signal, so that = ^^Wj.. Then the vector 
Sfc = GJyfc is the interference-free estimate of the original 
transmitted vector Sfc. 
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Colored Noise Caused by Interference 




Interference Channel 



Interference Channel 



Fig. 2. Receivers modeled by an interference channel may experience 
uncoordinated interference, modeled as colored noise, from part of the 
network not modeled as being a part of the same interference channel. 



B. Minimum Interference Plus Noise Leakage (INL) 

The subspace approach of |13|, p3) , outlined in Sec- 
tion 



III-A 



aims at aligning interference, which is capacity- 
optimal as the ratio of signal power to receiver noise power 
tends to infinity. If colored noise exists in any receiver, how- 
ever, the lA subspaces might be chosen to align with the noise 
to cancel it as well as the interference. Such colored noise may 
be due to an interference source outside of the coordinated 
portion of the network modeled as an interference channel, 
as shown in Figure |2] This interference is referred to as 
uncoordinated interference. We therefore focus on algorithms 
that take noise into account in their optimization. Note that 
these approaches have a "global" objective function limited to 
the users cooperating in interference channel, such as inside a 
single cluster in Figure |2] and thus assume the uncoordinated 
interferers of other clusters have fixed covariance over the 
optimization and transmission time. 



The objective of the subspace algorithm of Section III-A 



IS 

to minimize the total post-processing coordinated interference 
power, also known as interference leakage or interference 
power in (13], (14]. Thus, one intuitive solution is to minimize 
the total interference plus noise leakage, or INL. Mathemati- 
cally, this is represented with the global performance function 



K 

E 

fe=i 



E 




(13) 



where v^, is the received noise vector observed at receiver k. 
Expanding the expectation and exploiting the independence of 
the signal and noise vectors, the objective becomes 



(14) 



k=l i=l 



where = Ev^v^ is the covariance matrix of the noise at 
receiver k. The objective is then 



minimize JiNL({F£},{*fe}) 
subject to F^Ff = f^I,^e {!,..., i^} 



(15) 



The constraints on the precoders and receive subspaces are 
identical to those in Section |III-A| Further, since JTinl is 
rotation-invariant to each of the variables, the solutions lie 
on the Grassmann manifold and techniques derived for it can 
be used. Since ||A|||, = tr(AA*), JTinl can be rewritten as 

K K 

'^INL = E E (Hfc.£F,F,*H*_^ + Rfc) , (16) 

1=1 

which for fixed {F^} is minimized by p6) 



/ 



k mm 



K 



Rfe + ^ Hfc^^F^F^H^ £ 



(17) 



V 



1=1 

t^k 



For the precoders {F^}, it is sufficient to note that, for fixed 
minimizing JTinl with respect to {Ff} is equivalent 
to minimizing JTia with respect to {F^}, as is seen by 
comparing ( 14 1 and (j9]). Thus, the precoder solution is identical 
This solution effectively tries to align the coordinated 
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interference with the dominant directions of the noise (or 
uncoordinated interference) if the noise has significant energy. 
In particular, if the noise is highly correlated spatially with 
a rank-one covariance matrix, then R^, = (T^.aj,a^ and this 
algorithm will attempt to align the interference to if 
possible. Such noise, which may correspond to a single-stream 
uncoordinated interferer not part of the cooperating network, 
might then be mitigated, although full removal is unlikely. We 
can also prove the following quantitative conclusions. 

Proposition 1: If Rfe = crp V/c £ {!,..., X}, 
minimizing JTinl is equivalent to minimizing JTia- 
Proof: By definition, 

2 



then 



JiNL 



K 

E 

k=l 
K 

E 

k=l 



= JlA 



kit: I 



l^k 



k.ti: e 



+ tr(*^Rfe*fe) 



K 

E 

k=l 



alSk- 



(18) 



Since the summation in ( 18 1 is independent of any of the free 



variables, minimizing JTinl is equivalent to minimizing Jia 
when the noise is spatially white at each receiver ■ 

Proposition 2: As pfe — >■ oo or ||Rfc||i? — >■ for all k, 
jTiNL converges to J^ia- Thus, the subspace algorithm with 
noise consideration has the same SNR scaling as the pure 
interference alignment algorithm. 

Proof: Define Afc as the largest eigenvalue of Hermitian 
matrix Rfe. Then, 



JiNL 



JiA 



K 

E 

k=l 



tr (*^Rfe*fe) 



< Jia + AfcS'fe. 



(19) 



k=l 
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For any arbitrary {Rfc}, we define a sequence of functions 



1.^ 



AfcS'fe, 



(20) 



fc=i 



corresponding to a sequence of noise covariance matrices 
any e > 0, 



= Rfc/n, so that ||Rfc||i? as n oo. Then for 



< e 



(21) 



for all n > Y.k=i ^kSk/e- ■ 
From the proof, we also note that JTinl > JiA and 
minjTiNL — iff Rfc is singular for all k and the columns 
of the interference-aligning receiver matrices {^k} lie in the 
null spaces of their respective noise covariance matrices {Rfc}. 
The metric JTinl is, in fact, likely to have a positive global 
minimum unless the total number of streams is reduced below 
the degrees of freedom of the network, even if the noise is 
correlated, because the noise subspaces at different receivers 
will be not perfectly alignable almost surely. Adapting the 
number of streams in the network to improve finite-SNR 
performance is an interesting problem that is beyond the scope 
of this paper. 

In an idealized system with Gaussian signaling, colored 
noise may correspond to uncoordinated interference from out- 
side the network of interest. For instance, consider the scenario 
of a cellular network across a metropolitan area. The strategy 
for this network may be to coordinate three adjacent sectors 
to use interference alignment (via subspace optimization) to 
transmit to one mobile per sector in the downlink. For a regular 
lA solution, the uncoordinated interference arriving at each 
receiver from sectors outside the coordination area would be 
ignored or modeled as spatially white. The min-INL algorithm 
would be able to exploit the knowledge of this uncoordinated 
interference and account for it as necessary. 



The algorithms of Sections III-A and III-B aim to align the 
coordinated interference, which in turn maximizes capacity in 
a fully connected high-SNR network. We have seen in Propo- 
sition [T]that in finite-SNR environments, white Gaussian noise 
does not change the solutions of subspace methods. Although 
this section has presented an approach for networks with col- 
ored noise, algorithms with better throughput performance in 
finite-SNR regimes are desired, especially since most networks 
are not likely to be fully connected and thus may operate with 
low interference-to-noise ratio (INR), where subspace methods 
are not likely optimal even with colored noise considerations. 
To illustrate the problem of implementing subspace algorithms 
in a real network, consider the following argument. Suppose 
all interfering links {Hfc.^}, k £ have a path loss coefficient 
/3 whereas direct links have a path loss coefficient of 1. The 
subspace precoder design will then not depend at all on the 
value of (3 since the scalar multiplication does not change the 
direction of the signal. If the receivers use their interference 
suppression filters {Ufc} to cancel the interference, then the 
throughput of the system will be independent of f3. Thus, 
subspace algorithms treat weak and strong interferers equally, 
without exploiting the possible capacity gains available when 
interference is weakened. If no noise exists in the system. 



this is perfectly fine, since the receiver will still have an 
interference-free signal that it could decode perfectly. Re- 
alistically, however, a dynamic network would benefit from 
adapting its behavior to the relative interference energy. As 



shown numerically in Section IV the algorithms proposed in 



Sections III-C and III-D are more suited to such adaptation 



than the subspace method of Section |III-A[ 

C. Mean Squared Error Minimization 

A common metric for accounting for noise in linear 
receivers in wireless communication systems is the mean 
squared error For example, a zero-forcing linear MIMO 
receiver simply inverts the channel, and results in coloring 
and amplification of noise. An MMSE receiver balances the 
effects of noise with that of inverting the channel depending on 
the relative energy of each. This same concept can be applied 
to interference alignment, where the transmitter and receiver 
balance their wish to align the coordinated interference with 
the need for keeping the signal level well above the noise. 

Joint MMSE designs for MIMO channels have been stud- 
ied for years and have been applied to the point-to-point 
model ||T7l-|[l9| and the broadcast channel |^. The 

development for the interference channel is distinguished from 
previous work in that precoders and receivers need to be 
designed for multiple transmitters and receivers, rather than 
just the multiple transmitters or receivers as in the multi-user 
case, or a single transmitter and receiver in the point-to-point 
case. 



As opposed to objectives discussed in Sections III-A and 



|III-B| the MMSE directly designs the receive spatial filters 
{Gfc}. That is, the output of the product G^yfc is the estimate 
Sfc of Sfc, and the MMSE criterion minimizes the expected sum 
of the norms between each and Sfc for all fc, yielding the 
objective 



Jmse = VE||G^yfc-Sfc| 



fc=i 



(22) 



Substituting (|3]l for results in a global performance function 
of 



K 



k=l 



K 



Jmse = X! Hfc^fcFfcSfc+^ Hfc^fF^Sf+Vfc -Sfc 



l^k 



and an optimization objective of 

minimize Jmse ({F J, {Gfc}) 
subject to \\¥i\\],< pi,,^ K}. 



(23) 



(24) 



Expanding the expectation and simplifying, the optimization 
is equivalent to 



Hfe,fcFfcF^H; ^. 



minimize Yl,k=i (*^fc (-^fc 

-2M{tr(G^Hfc,fcFfc)} 
subject to < pt,le{l,...,K} 



(25) 



In general, MMSE solutions with an orthogonality con- 
straint are more difficult to derive. Thus, we relax the or- 
thogonality constraint to a total power inequality constraint 
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1 1 III' < pi, y£, and resort to a solution satisfying the Karush- 
Kuhn-Tucker (KKT) conditions as in previous joint MMSE 
solutions for different channel models [27 1. As shown in 
Appendix |A[ at each step the optimal precoders are 

-1 



\ k=l ) 



(26) 



where /i£ is the Lagrangian multiplier chosen to meet the 
power constraint. This may require a simple optimization 
(detailed in Appendix \A\ and has no known closed form. The 
optimal receivers are 




FTTi^Trr* 



(27) 



where no further optimization needs to be performed because 
there is no constraint on the receiver As the following propo- 
sition shows, this design can be viewed as a generalization of 
previous designs for the point-to-point case. 

Proposition 3: With tlk,e = for all £, k such that k ^ £, 



( [26] l and ( 27 1 are equivalent to an MMSE design for a point- 
to-point scenario. Further, as pk oo, the precoders and 
receivers diagonalize their respective information links. 

Proof: This is proven by substituting for each £, 
k £, and referring to previous point-to-point results p8) , 

m- ■ 

We also note that at high SNR and no coordinated inter-user 
interference, the MMSE algorithm will converge with one step, 
since any initialization precoder is a fixed point of the 
algorithm and will minimize the MSE. 

The MMSE design is unique among those discussed in this 
paper As discussed before, the MMSE receiver gives a direct 
estimate of s^, while the others require a conventional MIMO 
receiver after is applied. The MMSE receiver solution for 
fixed Fe,£ e {1, . . . ,K}, is simply the conventional MMSE 
MIMO receiver with colored noise. Further, the solution at 
each step is not in closed form, as an optimization needs to be 
done to meet the power constraint for the precoders. Lastly, the 
precoder solution is not orthogonal (or, conversely, a solution 
with orthogonal constraints is difficult to find). This algorithm 
may be difficult to implement because of these properties. 

Finally, the min-INL optimization is equivalent to an 
MMSE problem that compares the post-processing output to 
#^Hfc fcFfcS;; instead of simply s^. That is. 



K 



(28) 



fe=l 



Thus, the receiver is expected to remove only the effects 
of coordinated interference and white noise instead of having 
to correct for distortion created by the channel as well. This 
output must then be sent to a MIMO equalizer to remove inter- 
stream interference before symbol-by-symbol demodulation. 

D. Signal-to-Interference-Plus-Noise-Ratio Maximization 



The original subspace algorithm presented in Section III-A 
minimizes post-processing coordinated interference energy. 



The min-INL algorithm in Section [Tll-B | adds consideration for 
noise leakage as well, which can improve performance under 



colored noise. The MMSE solution in Section III-C indirectly 
accounts for signal power by attempting to force the received 
signal to look like the intended signal before precoding and 
transmission. It is clear, however, that a more desirable metric 
for maximizing the sum throughput would directly account 
for the post-processing signal-to-interference-plus-noise ratio 
(SINR). 

This section presents an algorithm for maximizing total 
SINR in the network. The optimization we use is not the 
only one that could be considered "maximum SINR", however, 
since total SINR for multiple nodes is not strictly defined 
in the literature. One may construct any number of global 
SINR metrics. Previous authors have considered the inter- 
stream interference for each transmit/receive pair and solved 
for the precoding and receiver matrices one column at a 



time 1 13 1, resulting in non-orthogonal precoders and receive 
spatial filters, as in the MMSE case. That approach, however, is 
not an alternating optimization of a global objective function, 
and its convergence is unproven. We therefore reformulate 
the problem into a maximization of the sum signal power 
across the network divided by the sum interference power, 
incorporating the inter-stream interference for each user The 
performance function becomes 



JsiNR — 



l^k=l l^n=l ''^l&fe 



e^k 



„(») 



Vfe 



(n) 

where is the nth column of matrix 



Eh.. 



(29) 



(30) 



is the pre-processing interference at receiver k from transmitter 



„(") _ V 

w — 1 



Hc(i") „w 
kMk ^k 



(31) 



is the pre-processing self-interference from streams w n at 
receiver k, and s)! is the nth entry of vector s^. Notice that 




Y^K Y^Si p(n) 



log (1 + Jsinr) 



(n) 

where P. is the post-processing signal energy of the 7ith 

(n) 

Stream at the ith receiver, / is the post-processing inter- 

(n) 

ference energy, and ' is the post-processing noise energy 



seen by the stream. The new objective (29 1 is the sum of 
signal power in the network divided by the sum coordinated 
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inter-user interference power and inter-stream interference 
power after processing. By maximizing this ratio the algorithm 
can design the precoders to either decrease post-processing 
interference (the denominator) or increase signal power (the 
numerator) to improve total network performance. 

The function JTsinr is a generalized Rayleigh quotient and 
can be solved using generalized eigen decomposition, and the 
optimization problem is 



maximize 



subject to 



||f{")||2 



in) 



},{gi"^} 

ynj. 



Se 



(32) 



For tractability we constrain each stream's precoder to have 
an norm equality constraint so that ||F^|||, — p(. For a larger 
objective function, and increased complexity, we could also 
introduce an inequality constraint on each column and vary 
the transmit power over the streams. As shown in Appendix [B] 
the solutions to the columns of the precoders are 

Se 



w = l 



.(to). 



H 



H 



fc=l ni= 



(n) 



(33) 



(n) 

where q\ is the sum of the terms in the denominator of (29 1 
that do not directly involve fl"'', and rl"'' is the sum of the 



p(n) 



terms in the numerator of ( 29 1 that do not directly involve 



The solutions to the columns of the receivers are 



(") _ I 



W — 1 



k,k 



K Si 

EE 

1=1 m=l 



Hf("l)f(™)*TT 



(n)p(n)* 



k.e 



(34) 



where i>max (A, B) is the generalized eigenvector correspond- 
ing to the largest generalized eigenvalue of the matrix pair 
(A,B). 



~(n) 
Ik 



and f 



in) 



are defined similarly as \33\ but with 

(") 



respect to ' instead of 

With all other variables fixed, the solutions in ( |33] l and (|34]) 
maximize the global SINR function ( [29] l, whereas the solutions 
in |13J give a suboptimal approximation to this solution. As 
shown in Section IV this does not imply that an iterative 



algorithm using the proposed solutions will converge to a 
larger objective than that of |13J. For any given channel 
realization and initialization, the two algorithms may give 
an identical result, or either may outperform the other The 



simulation results in Section IV suggest, however, that the two 
algorithms perform similarly on average. Since the proposed 
design requires more network knowledge than that of flj). 
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Fig. 3. Sum rate vs. initialization for each algorithm discussed in Section [ni| 
run on a single channel realization of the (2, 2, 3) MIMO IC at 10 and 40 dR 
Although the MMSE algorithm varies the most for this channel realization, 
this is not a aeneral trend. 



the latter is more attractive for implementation. If the extra 
network state knowledge is available, however, an intelligent 
design would be to run both algorithms and choose the design 
that works best for each channel realization, resulting in a 
sum throughput higher than either algorithm could produce 
individually. 

Note that the lA algorithm will minimize the left-hand term 
of the denominator in i7sinr, and the min-INL algorithm 
will minimize the entire denominator (minus inter-stream 
interference). Certainly, with no noise (or, more rigorously, as 
total signal to noise ratio goes to infinity), the two solutions 
are equivalent since maximizing the SINR will reduce to 
maximizing the SIR, which, as discussed before, lA does. This 
fact was proven in p4) . 

E. Convergence and Initialization 

This section analyzes some important details of the algo- 
rithms proposed in this paper In particular, the focus is on 
variable initialization, algorithm convergence, method of exe- 
cution, obtainment of channel state, and precoder constraints. 
We have found heuristically that arriving at a globally opti- 
mum point for the minimization algorithms (global optimality 
cannot be identified with the max SINR algorithm) is highly 
likely even when initializing the precoders to truncated identity 
matrices; the throughput, however, is the real objective we 
wish to optimize, and these algorithms only approximate that 
optimization. Thus, different initializations of an algorithm 
may result in drastically different throughputs, even if they 
result in the same final objective (or cost) function. For 
example, consider Figure [3] Each of the algorithms discussed 
in Section [Hi] was run on a fixed channel with 10 different 
random precoder initializations for the (2,2,3) MIMO IC at 
p = 40 dB and p = 10 dB. For p = 40 dB, the MMSE 
algorithm varied most between different initializations, but 
this is not indicative of the algorithms on the whole, just the 
behavior for this particular channel. It appears that finding 
"good" initializations is not difficult; experimentation has 
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shown that random initializations give as good of rates in these 
algorithms as any "intelligent" initialization tried. If possible, 
multiple runs of the algorithm should be made with different 
initializations for the best performance in terms of throughput, 
as shown in Figure |3] 



Each of the algorithms from Section III are guaranteed to 
converge because the objectives are bounded and at each step 
are moving monotonically in the direction of that bound. Con- 
vergence to a global optimum is not guaranteed except when 
the objective has certain convexity-like properties p2) that 
these algorithms are not proven to possess. Also, convergence 
of the cost function does not automatically imply convergence 
of the precoder designs, the analysis of which is beyond our 
scope. 

IV. Simulations 

This section presents simulations of the algorithms pre- 
sented in Section [Hi] to substantiate our claims and show that 
each of the algorithms can outperform the others in different 
regimes since none explicitly maximizes throughput. All of the 
simulations evaluate the expected sum rate with i.i.d. zero- 
mean unit-variance complex Gaussian coefficients for each 
channel, with the precoders for each realization calculated 
with perfect CSI and as if the realization was flat in time and 
frequency. More realistic channel scenarios are considered in 
our related work |28|. Transmitter k is assigned a deterministic 
transmit power pk and the link from transmitter £ to receiver 
k has a deterministic path loss coefficient ak.e- Whereas in 
preceding analysis a^^e was absorbed into H^.^, in this section 
we pull it out for exposition. We also define 7^,^ = oik,iPk.i 
to be the expected SNR at receiver k from transmitter I. Thus, 
the sum rate is 



' K 

E 



log 



I + R-fc ^afc,fcHfc_fcF/cF^H^ J. 



where 

R-fc = Rfe + ^ aA;^£Hfc_^F^F|H^ £ 

is the interference plus noise covariance. Precoders are ini- 
tialized randomly with orthonormal columns, as discussed in 
Section III-E and each algorithm is presented with identical 
initializations. Five random initializations are used for each 
channel realization, as motivated in Figure [5] and the initial- 
ization that maximizes ( 35 1 is kept while the others are thrown 
away. In each plot presented in this section, i?sum is computed 
via Monte Carlo simulations using 1000 independent channel 
realizations. Each iterative algorithm is run with 100 iterations 
each. 

Each algorithm from Section [III] is compared with a random 
precoding scenario where each precoder F^ is chosen as the 
left singular vectors of a random Gaussian matrix, to enforce 
an orthogonality constraint. That is, by Random Beamforming, 
we mean that 

F, = ,/^uf^\ (37) 

V 

where U^'^'^'' are the first S( columns of the left singular 



matrix of a random matrix with i.i.d. zero-mean unit-variance 
complex Gaussian coefficients. A greedy approach is also 
included to show the benefit of cooperation in the MIMO 
interference channel Q, |29|. In this design, each precoder 
Ffc, fc 7^ ^, is held fixed when designing F^. Then 



u.s.v: = R 



and 



a^,fcH£_fcFfcF^H|_j, 



'1/2 



(38) 
(39) 



The greedy algorithm is not guaranteed to converge since it is 
not optimizing a global function, but it requires less channel 
estimation. Finally, when the K — Z user interference channel 
is considered, the closed-form solution from |1| is also used 
for a baseline comparison. 

We first introduce colored noise into the interference chan- 
nel via an uncoordinated rank-one interferer in the network. 



as discussed in Section III-B Defining H^ ^; as the MIMO 



channel from the uncoordinated rank-one interferer to receiver 
k in the interference channel, then receiver k observes 



Yk 



K 

E 

1=1 



esi + y/ak,Eiik,EtESE +^k, (40) 



where since the uncoordinated interferer is rank-one, it is 
sending a single stream se precoded with vector f^;. Each 
receiver sees spatially white additive noise on top of the 
signal and interference (coordinated and uncoordinated). Pure 
interference alignment will ignore the uncoordinated interfer- 
ence, implicitly assuming it is spatially white. The rest of the 
algorithms will take the uncoordinated interferer into account 
but will not be able to fully suppress it without reducing the 
number of streams in the network since the uncoordinated 
interferer, which is scaling its power with the transmitters 
inside the network, is reducing the degrees of freedom of the 
network, making it interference-limited. 

Figure |4] illustrates results for the rank-one interferer sce- 
nario for each algorithm discussed in Section III with K ^ 3 
users, AI — N = 2 antennas at each node, and 5 = 1 
stream being transmitted between each transmit/receive pair 
for Pk — Pe — Pi Vfc and a = 1. That is, the transmit power is 
equal at all transmitters, including the uncoordinated interferer, 
and the path loss and fading statistics are identical on all links. 

The MMSE algorithm has higher degrees of freedom in 
this case because of its power inequality constraint on the 
precoders. This allows two transmitters to effectively shut off 
while the third has a degree of freedom and can cancel the 
external interferer with its extra receive antenna. This shows 
the flexibility of the MMSE design. Other than MMSE, the 
max-SINR algorithm outperforms the others in the power 
ranges considered. Note that, although on average the max 
SINR algorithm and the approximate max SINR algorithm 
have nearly identical performance, for any given channel re- 
alization they may have very different sum rates. I A performs 
the worst of all four iterative algorithms since it is neglecting 
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Fig. 4. Sum rate vs. = Pe = P for each algorithm discussed in Section [lIl] 
for the case where a rank-one uncoordinated interferer is introduced into the 
(2,2,3) networlc with 5 = 1 stream per user. The interferer's transmit power 
is scaled with the transmitters in the network so the degrees of freedom are 
reduced, and the network is interference-limited at high values of p. The 
MMSE algorithm is an exception because it has a power inequality constraint 
on its precoders and can thus allow two transmitters to turn off, giving the 
remaining transmitter one degree of freedom, so the sum capacity scales 
linearly with p. 

the uncoordinated interference. At high p, considering the 
colored noise in the algorithm objective results in a roughly 
20% increase in sum rate for this scenario. Note that the two 
best-performing algorithms, MMSE and max-SINR, do not 
have orthogonal precoders and thus may be more complex to 
implement in a real system with feedback requirements. With 
its orthogonal design and improved performance over lA, the 
min-INL algorithm is a good tradeoff between complexity and 
performance in this scenario. 

Next, we keep the same scenario but with fixed uncoordi- 
nated interference power, so that the degrees of freedom are 
not reduced. Figure |5] gives the results of this experiment. 
It shows that the uncoordinated interference, which is fixed 
at pe = dB, has little effect on the system, even at low 
Pk = p- The algorithms, except random beamforming, all 
scale at the same rate, and thus all exploit the maximum 
degrees of freedom in the network. For a fixed number of 
iterations, however, the MMSE algorithm does not scale, as 
it appears to require more iterations to converge than the 
others at high p. In particular, as shown in Figure |5] when 
the MMSE design is run with 500 iterations, its performance 
approaches that of the rest of the designs, while the other 
algorithms benefit very little from the increase in iterations. 
This is consistently seen in the rest of the simulations in 
this section. Analysis of this longer convergence is left to 
future work. Finally, we note that iterative lA outperforms the 
closed-form solution because multiple lA solutions exist, and 
iterative lA is better able to find the best one because of the 
multiple random initializations. If the closed-form algorithm 
is modified to explore multiple possible solutions, it would 
perform equally well in this case. 




p(dB) 



Fig. 5. Sum rate vs. pf^ = p for each algorithm discussed in Section [TlT| for 
the case where a rank-one uncoordinated interferer with fixed transmit power 
of = dB and pE = 20 dB is introduced into the (2,2,3) network with 
5 = 1 stream per user. The degrees of freedom in this network are the same 
as if the interferer did not exist, and each algorithm, with the exception of 
random beamforming, performs very similarly and exploits all the degrees of 
freedom in the network. 



30 r 
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p(dB) 

Fig. 6. Sum rate vs. = pE = p, with ^ = for fc > 1 and 
a-i E = li for E^ich algorithm discussed in Section llll| In this case, a rank- 
one uncoordinated interferer is sensed at only receiver 1 in the (2,2,3) MIMO 
interference channel with 5 = 1. The interferer's transmit power is scaled 
with the transmitters in the network so the degrees of freedom are reduced. 
The network is not interference-limited, however, since only one receiver sees 
the interference. 

Now we remove the uncoordinated interferer from all but 
one receiver in the network, and allow that uncoordinated 
interference power to scale with internal network transmit 
power, so that pk = pe = p, but ak^E = for fc > 1 
and ai^E — 1- Figure |6] shows the results. The maximum 
SINR and minimum INL algorithms suffer at high p relative 
to pure interference alignment and MMSE. This is because 
these algorithms see the large interferer at receiver fc = 1 as 
something to be overcome; these algorithms are effectively 
concerned with the average performance of the network. lA 
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Fig. 7. Sum rate vs. a = ai^ £, k ^ £, for the (2, 2, 3) MIMO interference 
channel with 5=1 stream per user. The iterative lA algorithm is not able 
to exploit reduced interference power, while all the other iterative algorithms 
can substantially improve throughput. The SNR on the data links is fixed at 
7fe,fe = 40 dB. 



is equally concerned about the average performance but only 
in terms of coordinated interference, whereas MMSE has 
flexibility to overcome the interference by reducing transmit 
power of receiver 1. To maximize sum rate in this case, it 
appears one should either ignore the external interference or 
include a power inequality constraint in the precoder design. 

We now turn to the case of no uncoordinated interfer- 
ence, considering only the conventional interference channel 
in isolation. In the first experiment, the transmit power is 
kept fixed but the path loss coefficient ak,e is varied on 
the interfering links (k ^ £) only. Figure |7] illustrates the 
results. The lA algorithm has constant throughput regardless 
of the interference path loss coefficient, but the other iterative 
algorithms are able to exploit the decrease in interference, 
converging to lA when the interference power is high. 

V. Conclusions and Future Work 

This paper has discussed the application and performance of 
iterative algorithms in the MIMO -user constant-coefficient 
interference channel under various operating regimes. The 
convergence and optimality of the algorithms has been dis- 
cussed, and similarities between all of them have been derived. 
If an iterative solution for the interference channel is ever prac- 
tical in a real system, it is unlikely that a direct interference 
alignment approach is desirable because of its suboptimality 
in environments where one or more links have little energy 
relative to the others. Instead, the max SINR or MMSE metrics 
are desirable in most environments because they flexibly adapt 
the solution between interference alignment (high interference 
power) and SVD precoding (no interference, fixed number of 
streams), and the MMSE solution in particular has a transmit 
power inequality constraint. These algorithms, however, have 
relatively high implementation complexity because of their 



nonorthogonality and lack of closed-form solutions at each 
step in general cases. In particular, the MMSE algorithm 
requires some optimization for meeting the power constraint, 
and the max-SINR algorithm requires more channel state 
knowledge at each iteration than the others. The min-INL 
algorithm is a good tradeoff between the three algorithms, 
since it has improved performance over lA in scenarios where 
there is uncoordinated interference or colored noise, but still 
has relatively low implementation complexity because of its 
simpler solutions and orthogonal precoders. Future work will 
focus on analyzing and reducing the overhead associated with 
solutions such as the ones presented in this paper. Although 
some studies have been carried out on the appUcation of in- 



terference alignment to a cellular network |30|-|32|, overhead 
and feedback analyses need to be performed to find out if the 
achievable gains are worth the effort. 

Appendix A 

Derivation of Mean Squared Error Minimization 

Proof: For completeness, we restate the optimization 
from (|24|, 



minimize Jmse ({Ff }, {Gfc}) 
subject to \\Fi\\j,<pe,£e{l,...,K}. (41) 

where Jmse is defined in ( p5| ). We use the Karush-Kuhn- 
Tucker conditions to solve the optimization at each step with 
all but one variable fixed. The Lagrangian of ( [24| is 

K 

2R {tr (G^Hfe^feFfe)} + ^ M£ (tr (F|F,) - pg) , (42) 

£=1 

where fig is the Lagrangian multiplier for the power constraint 
for precoder £. The KKT conditions are 



V/: = 

Ai£(tr(F,*F,)~ 1) = 0,V£ 

tr(F|F,) < pey£ 

fie > oyt 



(43) 
(44) 
(45) 
(46) 



For fixed {F^} and {fie}, {Gk} can be found by solving 
^G,C = Vg.Jmse = for fc e {!,..., if} and the 
KKT conditions will be automatically met since there are no 
constraints on {Gfe}. This yields 

Gk^ (j^Jlk^FiF^nie + Rk^ Hfe,fcFfe, (47) 



In solving for {F^}, we must ensure all of (43i-(46i are 
satisfied. To satisfy Vp^^C = 0, we must have 



Fe - ^fiel + ^ H.ljGkG*kFik,e^ 



H« cGf 



(48) 



If fie = satisfies ( [45] l, then all the KKT conditions are 
satisfied and the optimal Fe has been found for this step of 
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the alternating minimization. Otherwise, we must solve for nth column of the A;th receive combiner is then 



fie > such that HF^ 



pe to satisfy the KKT conditions. 



Although there is no known closed-form solution for pg in gj_' 
this case ||27|, HF^H^ is a monotonically decreasing function 



of He for lie > 0, so simple one-dimensional searches such as 
the bisection method can be done to solve for j^^}. ■ 



Appendix B 
Derivation of SINR Maximization 



(«) 



w — l 



K St 



e=i m=l 



p(n)H 



H 



k,k 



k,k 




(55) 



Proof: For completeness, we restate the optimization 
from (|32]i. 



maximize Jsinr ({f^"^}, {si"^} 
subject to 



f(")||2 



(49) 



where i7sinr is defined in (29 1. The optimization is performed 
on the columns of each precoder and spatial equalizer Gfe. 
In solving for the nth column of precoder Fe, we hold fixed 
every other precoder F^, k ^ £, and other columns f^^'"^ w ^ 
n of precoder F^, as well as every receive combining matrix 



Gfe, Vfc. The objective of (29 1 can be rewritten as 



Jsinr 



where 



EK v^Sfc f(n)* f pi(ri 
1=1 l^n=l '■£ \ye 



Se 



J") 



I + A 



B 



i 

(50) 



(n) 



B 



Si 

"-ejSe Se 

W — l 

K Sfc 

E\ ^ TT* „W„('")*-CT 
^k,iSk gfe ^k,l 

k=l m=l 



(")„(")*TT 
g£ til 



(51) 



(52) 



(53) 



(n) 

and ff^ is the remaining summation terms in the denominator 



of ( 29 1 that are independent of , contracted here for 



in) 

brevity. Similarly, r\ is the remaining summation terms 



in) 

in the numerator of (29i that are independent of . The 



function in ( 50 1 is the generalized Rayleigh quotient which is 



well known to be solved by the generalized eigen-vector of 
the numerator and denominator matrices. 



An) 




(n) 



B 



(») 



Ser'-rh 



(54) 



The derivation of the receive combiner columns follows the 
same structure as the precoders, but with a unit-norm con- 
straint on the columns for simplicity (this removes the Se 
multipliers in front of q^"'' and r^"^}. Here we define fL"'' 
to be the terms in the numerator of ( 29 1 independent of g^"'' 



in) 

and similarly for qj, ' in the denominator The solution for the 
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